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!♦  Introduction 


The  diffraction  of  guided  optical  waves  is  of  major  importance  in  integrated  optics 
owing  to  the  versatility  and  the  wide  use  of  integrated  periodic  structures.  Periodic 
interdigitated-electrode  devices  are  widely  used  in  integrated  optics  to  induce  phase  grat¬ 
ings.  These  phase  gratings  can  diffract  guided  modes  and  thus  can  function  as  switches*'*, 
modulators, deflectors,*  mode  converters,*"  **  strain  wave  transducers,*®  and  perform 
numerical  operations  for  optical  signal  processing**"**  and  optical  computing**"®*  appli¬ 
cations. 

Varioiis  methods  of  analyzing  the  guided-wave  diffraction  have  been  used.  The  most 
extensively  used  method  of  analysis  is  the  coupled-mode  approach  first  introduced  by 
Marcuse®*  and  Kogelnik®*  and  then  applied  to  waveguide  gratings.®^"**  Other  methods 
used  are  the  Floquet-Bloch  approach,*®"**  which  computes  the  exact  modes  in  the  per¬ 
turbed  region,  the  Green’s  function  method,**"**  first-order  perturbation  theory  applied 
to  the  total  field,**"*®  and  the  Rouard’s  method.**  **  In  almost  all  published  analyses 
mentioned  above,  the  interaction  of  the  optical  field  with  the  grating  is  collinear  which 
means  that  the  incident  wavevector  and  the  grating  vector  are  codirectional  or  contradi- 
rectional.  Most  of  these  analyses  apply  to  distributed  feedback  lasers,  distributed  Bragg 
reflector  lasers,  grating  couplers,  filters,  multiplexers  and  demultiplexers,  deflectors,  mode 
converters,  etc.  The  analyses  of  Refs.  33,  35, 51,  53,  and  54  can  also  treat  oblique  incidence 
on  pure  reflection  gratings  whose  grating  vectors  are  perpendicular  to  the  boundary  be¬ 
tween  the  waveguide  and  the  grating.  As  a  result,  only  one  transmitted  and  one  reflected 
beam  are  possible.  In  Ref.  36  an  approximate  ray  techmque  and  an  approximate  coupled- 
beam  technique  (taking  advantage  of  the  paraxial  beams)  are  used  to  treat  grating  vectors 
of  general  orientation  with  respect  to  the  incident  wavevector.  However,  only  the  Bragg 
diffracted  beam  is  considered  in  that  analysis.  Similarly  in  Refs.  47  and  48  the  general 
non-collinear  case  is  appraxicbately  treated  by  retaining  only  the  Bragg  diffracted  beam, 
using  Green’s  functions  or  first-order  perturbation  theory  respectively.  Kenan®*  has  ana- 
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lyzed  approximately  a  more  general  case  with  many  diffracted  orders  and  modes  for  both 
polarizations  using  a  coupled-mode/coupled-wave  approach.  In  all  the  above  analyses  the 
anisotropy  is  completely  ignored.  Even  though  the  coupled-mode  theory  has  been  extended 
to  anisotropic  media, its  application  to  non-collinear  interactions  is  cumbersome  due  to 
the  difficulties  in  the  computation  of  the  coupling  coefficients,  and  in  the  computation  of 
normal  modes  of  anisotropic  waveguides.  As  a  result,  very  little  work  has  been  done  on 
anisotropic  guided-wave  diffraction,  and  that  which  has  been  done  is  only  for  the  case  of 
collinear  interaction.^*  However,  anisotropic  materials  like  lithium  niobate  are  widely  used 
in  integrated  optics  and  consequently  the  anisotropic  properties  should  not  be  ignored. 

In  this  paper,  a  new  method  of  analyzing  the  diffraction  of  the  pure  guided  modes 
of  anisotropic  waveguides  by  inter  digitated-electrode  induced  phase  gratings  is  presented. 
The  analysis  includes:  (1)  the  anisotropic  properties  of  the  waveguide  and  of  the  grating, 
(2)  the  grating  vector  not  being  parallel  or  antiparallel  with  the  incident  mode  (non- 
collinear  diffraction  problem),  (3)  an  arbitrary  number  of  diffracted  orders  retained  in  the 
analysis,  (4)  the  conditions  for  mode  propagation,  (5)  the  conditions  for  efficient  mode 
diffraction  based  on  the  various  Bragg  conditions  of  the  plane  wave  components  of  the 
mode,  (6)  computation  of  the  mode  diffraction  efficiencies,  and  (7)  generalization  of  (l) 
through  (6)  for  the  treatment  of  the  guided-wave  diffraction  by  a  cascaded  stack  of  an¬ 
isotropic  gratings  induced  by  interdigitated-electrodes.  The  analysis  method  described  in 
Ref.  57  is  applied  for  the  computation  of  the  characteristics  of  the  induced  grating.  The 
diffraction  analysis  is  an  extension  of  the  rigorous  coupled-wave  approach**  based  on  the 
fact  that  a  pure  guided  mode  in  an  anisotropic  waveguide  can  be  decomposed  into  four 
homogeneous  plane  wave  components  (two  ordinary  and/or  two  extraordinary)  that  are 
not  phase-matched  at  the  grating  boundary  and  consequently  the  diffraction  problem  can 
be  decomposed  into  four  diffraction  subproblems.  The  analysis  is  restricted  to  uniaxial 
waveguides  with  one  principal  dielectric  axis  perpendicular  to  the  electrode  surface,  and 
to  electrode  periodicities  chosen  such  that  negligible  intermodal  coupling  occurs  (or  equiv- 
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alently  for  single-mode  waveguides).  The  general  approach  of  Ref.  59  has  been  adopted 
and  modified  in  order  to  compute  the  effective  mode  indices  and  the  fields  in  the  film 
region  of  the  waveguide  for  the  two  optic  axis  orientation  cases.  In  addition,  a  new  real 
transcendental  equation  for  the  computation  of  the  mode  propagation  constants  has  been 
derived  for  the  case  in  which  the  optic  axis  lies  in  the  electrode  plane.  The  conditions 
for  mode  propagation  are  discussed.  Those  are  classified  into  geometrical  conditions  and 
amplitude  and  phase  matching  conditions.  The  diffracted  mode  paursimeters  and  efficien¬ 
cies  are  also  calculated.  The  analysis  is  generalized  for  slanted  and  cascaded  electrode 
configurations.  It  has  been  found  that  the  most  efficient  configuration  occurs  when  the 
optic  axis  is  perpendicular  to  the  electrode  surface.  For  the  case  of  the  optic  axis  lying 
in  the  electrode  plane,  the  optic  axis  orientations  parallel  or  perpendicular  to  the  grating 
vector  are  treated.  Finally,  some  example  devices  are  analyzed.  These  include  a  Givens 
rotation  device,  a  herringbone  grating  multiplier,  and  a  hybrid  mode  diffraction  example. 
In  addition,  a  comparison  with  experimental  data  is  included  and  an  efficient  configuration 
is  suggested. 

The  geometry  of  the  problem  is  shown  in  Fig.  la.  In  Fig.  lb  a  top-view  of  the 
same  configuration  is  shown.  In  both  figures  (z,  y,  z)  is  the  coordinate  axis  system  that 
corresponds  to  the  three-dimensional  structure,  while  (z, ,  y, ,  z» )  is  the  coordinate  axis 
system  of  the  slab  uniaxial  anisotropic  waveguide  which  is  rotated  by  an  angle  S  about  the 
z  axis  of  the  (z,  y,  z)  coordinate  system.  The  analysis  is  based  on  the  given  conditions: 

(1)  The  optic  axis  is  restricted  to  He  in  the  electrode  plane  or  perpendicular  to  that  plane. 

(2)  The  intermodal  coupling  is  negligible.  This  is  automatically  true  if  these  are  single- 
mode  waveguides.  (3)  The  evanescent-field  diffraction  has  the  same  characteristics  as  the 
corresponding  homogeneous  plane  wave  diffraction.  (4)  The  fringing  effects  at  the  grating 
boundaries  are  neglected. 


n.  Guided- Waves  in  Uniaxial  Anisotropic  Slab  Waveguides 


The  modes  that  can  propagate  in  uniaxial  anisotropic  dielectric  slab  waveguides  are 
discussed  in  this  section.  The  optic  axis  orientations  in  the  film  and  the  substrate  regions 
are  the  same,  and  sje  restricted  by  the  condition  (l)  given  above.  The  cover  is  taken 
to  be  isotropic  (corresponds  to  the  buffer  layer  that  is  used  for  integrated  electro-optic 
applications  and  it  usually  consists  of  silicon  dioxide).  The  geometry  of  the  waveguiding 
system  is  also  shown  in  Fig.  la. 


The  tensor  relative  permittivities  of  the  film  and  the  substrate  are  £/ ,  and  f,  respec¬ 
tively  (expressed  in  the  (x, ,  y,,  x.)  coordinate  axis  system),  while  the  scalar  relative 
permittivity  of  the  cover  is  e*.  The  thickness  of  the  waveguide  (film  thickness)  is  d  and 
the  coordinate  axis  system  (x. ,  y. ,  x« )  is  related  to  the  (x,  y,  x)  coordinate  system 
by  a  simple  translation  and  rotation  transformation.  Plane  wave  solutions  of  the  form 
cxpl-i^o(*^««.*«  +  0'^w)]  used  where  i  =  c,f,s  correspond  to  the  cover,  film,  and 
substrate  regions,  respectively.  The  Jfc,*,-  and  are  the  normalized  wavevector  compo¬ 
nents  (0'  is  usually  known  as  the  normalized  propagation  constant  or  effective  index)  and 
ko  =  2n/Xo  where  Ao  is  the  freespace  wavelength.  Using  the  analysis  of  Refs.  59  and  60  the 
complex  tangential  electric  and  magnetic  field  components,  E,^i,  H,,,,  and  H,,, 


cam  be  written  in  the  form 


(1) 


where  JF,  =  is  a  4  x  1  matrix  (vector),  >4^  is  a  4  x  4  matrix 

given  in  Appendix  A,  i  =  e,/,s,  and  ijo  =  (#*o/co)^^*.  »  the  characteristic  impedance  of 
freespace,  where  Co»  *ind  fio  are  the  permittivity  amd  the  permeability  of  the  freespace. 
The  solution  of  Eq.  (1)  is 


Ei  =  [Wi  exp(-yfcoi«*.)l(^<  exp(-j7:o^'*.)].  (2) 


where  is  a  4  x  4  matrix  containing  the  eigenvectors  of  matrix  A^,  Ai  is  a  diagonal 
4x4  matrix  containing  the  eigenvalues  of  A,  ,  and  is  a  4  x  1  matrix  containing  the 
unknown  constants  of  the  first-order  differential  equations  [Eq.  (1)].  The  elements 
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(p,g  =  x,i/,z  and  i  =  c, /,«)  correspond  to  the  permittivity  elements  of  each  region  and 
are  expressed  in  the  (x, ,  y, ,  x, )  coordinate  axis  system.  The  dispersion  equation  for 
the  normalized  propagation  constant  /?'  can  be  found  using  the  boundary  conditions  of 
the  problem.  These  are  described  in  Appendix  A.  It  is  straightforw’ard  to  show  (using 
the  expressions  of  Appendix  A)  that  for  nonzero  solutions  of  the  tangential  fields,  the 
normalized  propagation  constant  p’  must  satisfied  the  equation 


P«(/?')P44(^')-P54(^')Po(^')=0, 


(3) 


where  p.y(/?')  is  the  tj-th  element  of  P(P'),  which  is  defined  in  Apppendix  A.  In  addition, 
for  evanescent  plane  wave  solutions  in  the  cover  and  substrate  regions,  and  homogeneous 
plane  wave  solutions  in  the  film  region,  the  normalized  propagation  constant  /?'  should 
satisfy  the  following  inequalities 


>  max{nl^,  n?>,, 

P 


3  3 

^MMw  $^stw  0  ^m*w0 


3  3 


(4) 

(5) 


where  no*,  n©/,  no,  are  the  ordinary  refractive  indices  of  the  cover,  the  film,  and  the 
substrate  respectively,  and  nf,/,  ng,  are  the  principal  extraordinary  indices  of  the  film 
and  the  substrate  respectively.  Inequality  (4)  guarantees  that  the  plane  waves  in  the  cover 
and  the  substrate  regions  are  evanescent  (both  ordinary  and  extraordinary  plane  waves) 
and  inequality  (5)  guarantees  that  both  ordinary  and  extraordinary  plane  waves  in  the 
film  are  homogeneous.  Inequalities  (4)  and  (5)  satisfy  also  the  power  flow  condition  since 
no  power  should  escape  normally  from  the  film-substrate  and  film-cover  interfaces.  The 
power  flow  condition  for  the  two  optic  axis  orientation  cases  that  are  considered  in  this 
chapter,  are  satisfied  because  the  Poynting  vector  components  normal  to  the  interfaces 
and  the  corresponding  wavevector  components  always  have  the  same  sign.  That  is,  > 
0  <=>•  >  0  and  5...  <  0  <=>>  <  0,  for  t  =  c,/,«  (5..<  is  the  component  of  the 

Poynting  vector  normal  to  the  interface).  This  condition  is  important  since  the  power  flow 
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in  the  cover  and  substrate  regions  should  be  zero  normal  to  the  interfaces  for  guidance  of 
the  light. 


If  inequalities  (4)  and  (5)  are  satisfied  then  all  plane  wave  components  in  the  film 
region  are  homogeneous.  These  homogeneous  plane  waves  for  which  0  satisfies  Eq.  (3) 
constitute,  along  with  the  corresponding  fields  in  the  cover  and  the  substrate,  a  subset  of 
the  possible  guided  modes  in  the  uniaxial  anisotropic  waveguide.®®-  In  this  analysis, 
only  these  kinds  of  guided  modes  are  considered.  This  restriction  limits  the  angle  of 
incidence  S  when  the  optic  axis  lies  in  the  electrode  plane.  Consequently,  (assuming  a 
negative  birefringent  material)  the  following  inequality  should  hold 


no,  < 


(6) 


f.../ (%.../ (^) 

where  all  the  permittivity  elements  are  expressed  in  the  waveguide  system  (x, ,  ) 

and  consequently  depend  on  the  angle  of  incidence  6.  In  the  case  of  positive  birefringent 
materials  the  inequality  for  the  angle  of  incidence  6  b  included  in  the  following  inequality 


3  ^  _ 


(7) 


where  in  both  cases  it  is  assumed  that  the  substrate  ordinary  and  principal  extraordinary 
refractives  indices  are  greater  than  no*.  In  the  following  subsections  the  two  optir  axb 
orientations  under  consideration  are  treated  separately. 


A.  Optic  Axis  Perpendicular  to  the  Electrode  Surface 

In  thb  case  the  optic  axb  b  oriented  along  the  z  axb  (Fig.  la).  The  modes 
are  decoupled  since  the  ordinary  and  the  extraordinary  waves  can  be  analyzed  sepa¬ 
rately.  Consequently,  TE  and  TM  modes  can  be  dbtingubhed.  Equation  (3)  becomes 
P33(^')P44(^')  =  0-  The  condition  Ps»(/9)  =  0  corresponds  to  the  guidance  condi¬ 
tion  for  the  TE  modes  and  b  given  by  2kokg^fod  ~  24>/co  ~  —  2m7r,  where 
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^••/O  —  (*^0/  ~  ^  *)*^*»  ^leO  —  ^  {|^«»e|/^*»/0  }»  ^/»0  ~  tail  ^  { |^«»  »£>  )/^»i»/0 }» 

ik„.  =  -y(^'’  -  n*  J'/*,  and  *,..0  =  “>{/?'’  -  The  integer  m  takes  the  values 

0,  1,  2,...,  and  is  characteristic  of  the  TE„  mode.  The  field  expressions  can  be  found 
using  Eq.  (2).  In  this  work  only  the  expressions  for  the  fields  in  the  film  region  aire  going 
to  be  presented  since  only  these  are  needed  in  the  remaining  analysis.  Those  expressions 
are  given  in  Appendix  A  [Eqs.  (A5)-(A8)]. 

Due  to  the  guidance  condition,  Eq.  (3),  one  parameter  of  the  problem  can  be  chosen 
arbitrarily  (depending  on  the  initial  conditions  of  the  problem  as  for  example  given  by 
the  input  power).  This  parameter  is  chosen  to  be  So-  (Eq.  (A4)].  Consequently,  the 
coefficients  Fo^ ,  and  Fo  -  (complex  in  general)  are  not  independent,  and  their  ratio  differs 
for  each  TE  mode.  These  coefficients  can  be  written  as  Fo+  =  A©  [€xp(-y^/,o)]5o- , 
Fo-  =  Ao (exp (+ji'0y.o)l‘S'o-i  where  Ao  =  (l/2)(l  +  |A:,.,o has 
been  previously  defined. 

The  condition  p*i{P')  =0  corresponds  to  the  guidance  condition  for  the  TM  modes 
and  is  given  by  2kok,^fgd  -  2^jtB  —  2d>f,E  =  2m?r,  where  k,,fB  —  {^oi ~ 
4>fcB  =  tan"*{n’ 4>f,s  -  tan"*{n^,lfc,,,«l/n* 
and  =  -j(no, /n£,)(^'’  -  n|,)*/^.  The  integer  m  takes  the  values  0,  1,  2,...,  and 

is  characteristic  of  the  TM^,  mode.  The  field  expressions  can  be  found  using  Eq.  (2). 
Similarly  to  the  TE  mode  case  the  fields  only  in  the  film  region  are  needed  and  are  given 
in  Appendix  A  |Eqs.(A9)-(A12)]. 

For  the  TM  mode  case  the  free  parameter  is  chosen  to  be  S'*-  [Eq.  (A4)).  Con¬ 
sequently,  the  coefficients  Fg^,  wid  (complex  in  general)  are  not  independent, 

and  their  ratio  differs  for  each  TM  mode.  These  coefficients  can  be  written  as, 
=  A*(exp(-j>,.*)]5£_,  Fg.  =  Aelexp(-I-;^„*)]S*-,  where  A*  =  (1/2)(1  + 
no/l*...r  previously  defined. 
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B.  Optic  Axis  Lying  in  the  Plane  of  the  Electrodes 


In  this  subsection  the  optic  axis  lying  in  the  xy  plane  is  treated.  The  general  theory 
of  Refs.  59  and  60  is  again  being  used.  The  guidance  condition  for  the  hybrid  modes  in 
this  case  is  derived  for  the  first  time  In  the  form  of  an  analytical  transcendental  equation 
similar  to  the  TE  and  TM  modes.  This  condition  is  given  in  general  by  Eq.  (3)  which  for 
this  optic  axis  orientation  becomes 

Go  Bo  Foo  E  +ho)sin(^o£  ~^oo)+ 

Ge  He  Fe  E  Fe  o  +h£)sin(^££  ~  ^eo)~^ 

Go  He  FooFeb  [siii(oio  ois  "h  ^oo  ■+■  s  +  Jo  ~  )~ 

sin(Q:o  —  ot£  +  4>oo  ~~  4*bs  9o  h®  )J— 

Go  He  Fo  B  Fe  ^os  +Jo  ~^b)~ 

sin(oio  ~  ctE  +^os  ~^BO  +  Jo  +h£)]+ 

Ge  Ho  Fo e  Fe  o  (sin(Q:o  +  ctg  +  4>eo  "h  ffs  ~  ^o)~ 

8m(a£  —  ao  +  —  ^on  +  Js  +  ^o)]— 

Ge  Ho  FooFe  e  [sin(Q;o  "h  "I"  ^oo  "i"  ^eb  9e  ~  )~ 

bui(q;£  —  Oco  +  (^B  B  ~  o  "h  j£  +  ho )]  =  0  (8) 

where  Fooy  Foei  Feo^  ■^se>  Go,  Ge,  Hq,  He,  <^0,  ^b,  9o,  9b,  ^o,  h^,  4>oo,  4>oe, 
4>bo,  ^<1  4>es  defined  in  the  Appendix  6.  All  the  above  parameters  are  functions  of 
the  normalized  propagation  constant,  fi'  {which  is  the  unknown  of  Eq.  (8)),  the  ordinary 
and  principal  extraordinary  refractive  indices  of  the  cover,  the  film  and  the  substrate 
regions,  the  thickness  of  the  waveguide,  'f,  the  orientation  of  the  optic  axis  and  the  angle 
of  incidence.  The  fields  in  the  film  region  can  be  found  in  a  similar  manner  to  that 
previously  presented  and  are  given  in  Appendix  B  [Eqs.  (B10)-(B15)].  Similarly  to  the 
decoupled  TE  and  TM  cases,  the  coefficients  Fo  + ,  Fo-,  Fs  + ,  and  Fe  -  that  appear 
in  the  fit  Id  expressions  are  not  independent.  Again  there  is  the  poisiollity  of  one  free 
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parameter,  which  as  in  the  previous  cases  is  chosen  to  be  So  -  •  Consequently  the  following 
expressions  can  be  written:  Fo+  =  Bo  (exp(-ytfo)]5o-,  Fo-  =  ■Bo[exp(+jiJo)|So-, 
Fe^  =  B£[exp(-jT>£)]So-,  F£_  =  B£ [exp(+jft?£)]So where  Bo,  iJo,  and  !?£  are 
given  in  Appendix  B. 


III.  Electro-opticallv  Induced  Phase  Grating 

When  a  voltage  is  applied  to  the  interdigitated  electrodes  (usually  a  dc  voltage  or  a 
quasi-static  voltage),  a  phase  grating  is  formed  inside  the  waveguide  under  the  electrodes. 
This  is  due  to  the  linear  electro-optic  effect  (Pockels  effect).  The  induced  grating  is  periodic 
along  the  z  direction  (Fig.  la)  and  decays  exponentially  along  the  z  direction  since  the 
inducing  electric  field  has  the  same  properties.  Consequently,  using  the  expressions  of 
Ref.  57,  the  optical  relative  permittivity  tensor  elements  in  the  waveguiding  region  (film) 
are 

(^,  ^)  —  "I"  •^m  Pm  (a)^m  (®),  (9) 

m,  odd 

where  u,v  =  z,y,  z  and  e»,o  is  given  in  Ref.  57,  is  a  constant  factor  that  depends  on  the 
electrode  geometry,  the  material  parameters,  the  applied  voltage,  and  the  spatial  harmonic 
number  m.  The  constant  is  given  in  Ref.  57.  In  addition,  ym(a)  =  exp[-7rma4(z  - 
t  —  de)/A,],  where  2t  is  the  electrode  thickness,  de  b  the  buffer  layer  (cover)  thickness, 
and  a4  is  also  defined  in  Ref.  57.  Finally,  ^m(z)  =  Bin(»rmz/A,)  where  A,  is  the  electrode 
spacing.  The  above  relative  optical  permittivity  tensor  elements  do  not  vary  along  the  ±y 
direction  and  are  expressed  in  the  (z,  y,  z)  coordinate  axis  system. 

In  order  to  account  for  the  variation  along  the  z  direction,  two  approaches  can  be 
used;  (a)  averaging  of  the  electric  field  (electrostatic  or  quasi-static  field)  with  respect  to  z 
coordinate,  without  any  optical  field  weighting,  and  (b)  averaging  of  the  electric  field  with 
respect  to  the  z  coordinate  using  optical  field  weighting.  The  inclusion  of  an  optical  field 
weighting  takes  into  account  the  non-uniformity  of  the  optical  field  in  the  film  region  of 
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the  waveguide.  Thebe  averaging  procedures  define  a  new  quantity, 

=  /  Pm  {2)dz,  (10) 

»  »0 

which  is  independent  of  For  the  case  of  no  optical  field  weighting, 

t"  =  /'■  W|£„.(*)l’<!»/  /'■  (11) 

Jmq  Jbo 

where  *0  =  *  +  <fe.  =  t  +  d*  +  d,  d  is  the  film  thickness,  and  \Eopt{z)\^  =  \E,^f  1’  + 
\E^„f\^  +  is  the  optical  field  magnitude  u  =  x,y,z  are  given  in  Eq.  (A5) 

or  Eqs.  (A9)  and  (A  10)  when  the  optic  axis  is  perpendicular  to  the  electrode  surface  and 
by  Eqs.  (bl0)-(bl2)  when  the  optic  axis  is  in  the  electrode  surface].  Since  the  optical  fields 
are  expressed  in  the  waveguide  coordinate  system  (z. ,  y. ,  2^),  a  suitable  coordinate  axis 
system  transformation  is  necessary  to  express  all  optical  fields  of  the  incident  mode  in  the 
(z,  y,  z)  coordinate  axis  system.  Substituting  Eqs.  (10)  or  (11)  into  Eq.  (9)  the  following  z 
independent  approximate  expressions  for  the  relative  optical  permittivity  tensor  elements 
are  found; 

(®)  —  ^mPm^m(2^)»  (l2) 

•n,  odd 

where  c**  (z)  is  the  uv-th  z  independent  element  of  the  relative  optical  permittivity  ten¬ 
sor.  Knowledge  of  the  relative  optical  permittivity  tensor  is  equivalent  to  knowledge  of  the 
electro-optically  induced  phase  grating  under  the  inter  digitated  electrodes  inside  the  wave¬ 
guide.  This  relative  optical  permittivity  tensor  is  necessary  to  carry  out  the  diffraction 
analysis. 


rv.  Guided  Mode  Decomposition 

Using  the  expression  of  Eqs.  (A5)-(A8),  or  Eqs.  (A9)~(A12),  or  Eqs.  (B10)-(B15), 
the  total  incident  optical  electric  and  magnetic  fields  are 

Eine  =  Oo+  (j*)  +  Eo-*io~[x)  +  Eg  +  as+i*)  +  Ea-O*-  (*)]6(z,y),  (13) 

Hinc  =  IH0  +  O0+  (*)  +  Ho-fio~{x)  +  +  a£  +  (2)  +  _  Oa_  (2)]6(z,y),  (14) 
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where  a,±(z)  =  exp(±j7:oA:.«/,*),  9  =  O,  JE?,  b[x,y)  =  exp(-jfco/?'(sin5  z  -  cos^  y)l, 
and  S  is  the  angle  of  incidence  (Fig.  la).  The  terms  (9  = 

O,  E)^  correspond  to  the  vector  complex  constants  of  the  incident  field  components.  The 
components  of  the  Eqs.  (A5)-(A8),  (A9)-(A12),  and  (B10)-(B15)  have  been  transformed 
into  the  (z,  y,  z)  coordinate  axis  system.  Using  Eqs.  (13)  and  (14)  a  guided  mode  incident 
on  the  grating  can  be  decomposed  into  four  plane  wave  components,  two  ordinary  0+  and 
O— ,  and  two  extraordinary  E+,  and  E—.  This  situation  is  depicted  in  Fig.  2. 

The  four  plane  wave  decomposition  corresponds  to  the  general  hybrid  mode  case. 
This  applies  when  the  optic  axis  is  in  the  plane  of  the  electrodes  (zy  plane).  For  the 
decoupled  cases,  when  the  optic  axis  is  perpendicular  to  the  electrode  surface  (along  2 
axis)  only  two  out  of  the  four  plane  waves  exist  and  the  mode  is  either  TE  or  TM.  For 
the  TE  decoupled  mode  case  only  the  0+  and  O—  plane  waves  exist,  and  consequently, 
Ee^  =  Ee-  =  He =  He-  =  0.  For  the  TM  decoupled  mode  case  only  the  E+  and  the 
E~  plane  waves  exist  and  consequently  Eo+  ~  Eq-  —  Ho*  —  Ho-  =0.  In  this  analysis 
it  is  required  as  mentioned  before  to  have  only  homogeneous  plane  wave  components  in 
the  film  region.  As  a  result,  inequality  (5)  should  hold.  As  described  in  Refs.  59  and 
60  it  is  possible  for  a  guided  mode  to  consist  of  both  homogeneous  and  inhomogeneous 
plane  wave  components  in  the  film  region.  Thus,  only  a  subset  of  the  guided  modes  in 
uniaxial  anisotropic  waveguides,  the  subset  that  contains  the  modes  decomposable  into 
four  homogeneous  plane  wave  components  in  the  film  region,  will  be  analyzed  using  the 
approach  described  in  this  section. 


V.  Decomposition  of  the  Diffraction  Problem  into  Four  Subproblems 

A  guided  mode  that  is  composed  (in  the  film  region)  of  four  homogeneous  plane  waves 
can  be  expressed  by  Eqs.  (13)  and  (14),  where  Ic.v/oi  KmfSt  >Ad  /?'  are  real  numbers. 
These  four  plane  wave  components  are  not  phase  matched  on  the  boundary  of  the  grating 
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[xz  plane  in  Fig.  la)  since  in  general  k,^fo  ^  incident  electric  and  magnetic 

fields  can  be  written  as 

Einc  exp{-yA:o  *«  •  *0 »  (15) 

« 

S,»c  •  f),  (16) 

« 

where  q  =  0+,  O— ,  E+,  E—,  and  k,  =  /?'sind  x  —  0'cos6  y  ±  k,^f^z.  The  +  sign 
corresponds  to  the  0+  and  waves  and  the  —  sign  corresponds  to  the  O—  and  E— 
waves.  Each  component  of  Eqs.  (15)  and  (16)  satisfies  the  Maxwell’s  equations.  Inside  the 
grating  region  the  total  field  is  expanded  into  space-harmonic  components  by  extending 
the  expressions  of  Ref.  58.  Those  expressions  are  extended,  however,  for  multiple  incident 
waves.  The  resulting  fields  in  the  grating  region  are 

=  X]  +  SI  (y)5]5.  (x,y)a,  (a),  (17) 

<  • 

S2  =  +  (y)y-i-^^*  (y)2l^(*,y)a.(«)»  (is) 

«  . 

where  q  =  O-H,  0-,  E-\-,  E-,  6i(x,y)  ^  exp{— yi:o((^' sind  -  iiif,)x  —  iK,y]},  and 
t  designates  the  t-th  diffracted  order  corresponding  to  the  9-th  incident  homogeneous 
plane  wave,  K  =  K,x  +  K^y,  is  the  grating  vector,  and  a,(a)’s  have  been  previoxisly 
defined.  The  E^  and  fields  should  satisfy  the  Maxwell’s  equations  in  the  grating  region. 
Substitution  of  Eqs.  (17)  and  (18)  into  the  Maxwell’s  equations  leads  to  four  independent 
sets  of  equations,  one  for  every  incident  plane  wave  component  (two  sets  of  equations  for 
the  ordinary  waves  and  two  sets  of  equations  for  the  extraordinary  waves).  The  above 
decomposition  of  the  Maxwell’s  equations  into  four  independent  sets  of  equations,  b  a 
result  of  the  linear  independence  of  the  exponential  terms  a,(s). 

El  =  Einc  +  exp(-y*ofcJ<  •  ^  txp{-jkok[\  •  f),  (19) 

*  i 

7?  exp(-y*of;<  ♦  (r  +  sy)]  +  7*’  exp[-jkok'^^  •  (r -h  sy)],  (20) 

f  i 
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where  Eine  is  given  by  Eq.  (15),  q  =  0+,  O— ,  E+,  E— ,  E*,  2^,  2^*,  are  the 

complex  diffracted  amplitudes  of  the  t-order  that  correspond  to  the  q  incident  plane  wave, 
and  /t’.,  k'^.  (p  =  1,3)  are  given  by  expressions  similar  to  the  ones  described  in  Ref.  58 
for  each  incident  plane  wave  q.  Using  exactly  the  same  procedure  described  in  Ref.  58 
for  each  incident  plane  wave  component  q,  four  independent  sets  of  boundary  conditions 
can  straightforwardly  be  derived.  Consequently,  the  three-dimensional  anisotropic  grating 
diffraction  analysis  of  Ref.  58  can  be  applied  repeatedly  for  each  incident  homogeneous 
plane  wave  component  of  the  incident  mode.  In  this  analysis,  the  condition  (3)  described 
in  section  I  is  being  used  since  the  boundary  conditions  for  the  film-cover  and  film-substrate 
interfaces  are  not  explicitly  used  inside  the  grating  region.  A  similar  condition  was  used 
by  Marc\ise’‘  in  the  coupled-mode  analysis  of  isotropic  corrugated  waveguides. 

Summarizing,  the  grating  diffraction  problem  of  an  incident  guided  mode  in  the  uni¬ 
axial  slab  waveguide  can  be  decomposed  into  four  plane  wave  grating  diffraction  subprob- 
lems.  Each  subproblem  corresponds  to  a  plane  wave  component  of  the  incident  mode.  The 
procedure  described  is  shown  in  Fig.  3. 


VI.  Mode  Propagation  Conditions 

A  pure  guided  mode  in  a  uniaxial  dielectric  slab  waveguide  can  be  decomposed  into 
four  homogeneous  plane  wave  components.  Thus,  in  general,  one  guided  mode  consists  of 
both  ordinary  and  extraordinary  waves.  Consequently,  all  types  of  diffraction  described 
in  Ref.  58  are  possible.  The  00,  OE,  EO,  and  EE  types  diffraction  are  due  to  the 
grating  anisotropy  and  to  the  hybrid  character  of  the  incident  mode.  Since  the  forward- 
and  backward-diffracted  wavevectors  satisfy  different  dispersion  equations  (ordinary  or 
extraordinary),  and  the  incident  wavevectors  are  not  phase  matched  along  the  z  axis,  the 
resulting  diffraction  angles,  projected  onto  the  xy  plane,  are  not  equal.  Thb  happens  since 
every  diffracted  mode  consists  of  four  ordinary  and  four  extraordinary  plane  waves  not 
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necessarily  lying  in  the  same  vertical  plane.  This  situation  is  depicted  in  Fig.  lb  for  two 
forward-diffracted  orders.  The  eight  plane  wave  components  result  from  the  diffraction  of 
each  incident  plane  wave  (mode  component)  into  an  ordinary  and  an  extraordinary  plane 
wave.  Since  the  incident  mode  contains  two  ordinary  and  two  extraordinary  waves,  eight 
plane  waves  result  from  the  PQ  types  of  diffraction  where  Q  =  E.  Due  to  the 
multiple  diffraction  types  there  are  corresponding  multiple  Bragg  conditions.** 

In  order  for  the  forward-  and  backw2Lrd-diffracted  modes  to  propagate  the  following 
conditions  must  be  satisfied:  (l)  At  least  four  of  the  eight  plane  wave  components  of  every 
diffracted  order  should  lie  on  the  same  vertical  plane  in  order  to  constitute  the  four  plane 
wave  components  of  the  mode  (two  of  the  plane  waves  should  be  ordinary,  O+t  0-,  and 
two  extraordinary,  E+,  E—).  For  the  decoupled  cases,  when  the  optic  axis  is  along  the  z 
axis,  only  two  plane  waves  are  necessary  to  satisfy  the  above  condition  (only  two  ordinary, 
0+,  0-,  in  the  case  of  TE  modes  and  only  two  extraordinary,  jP-f ,  jP—,  in  the  case  of 
TM  modes).  (2)  The  propagation  constant  of  a  diffracted  mode  should  be  the  same  as  the 
propagation  constant  of  the  incident  mode  along  the  diffraction  direction  that  corresponds 
to  that  particular  diffracted  order  (since  in  the  uniaxial  dielectric  slab  waveguides  the 
mode  parameters  depend  on  the  direction  of  the  propagation  with  respect  to  the  optic  axis 
orientation).  This  condition  is  necessary  to  guarantee  that  the  same  type  of  mode  will 
propagate  in  the  uniaxial  slab  waveguide  after  the  diffraction  by  the  induced  phase  grating 
(since  either  the  intermodal  coupling  is  neglected  or  the  waveguide  supports  only  a  single 
mode).  (3)  For  efficient  diffraction  the  same  Bragg  conditions  should  be  satisfied  by  all  the 
incident  plane  wave  components  of  the  mode.  (4)  Finally,  the  four  plane  wave  components 
of  the  mode  should  satisfy  a  certain  amplitude  and  phase  condition  which  is  characteristic 
of  the  incident  mode.  In  the  case  of  hybrid  modes  this  condition  will  catise  radiation  losses 
since  in  general  the  ordinary  and  the  extraordinary  components  of  the  mode  will  possess 
different  diffraction  characteristics. 

Conditions  (1),  (2),  and  (3)  restrict  the  direction  of  the  optic  axis  with  respect  to 
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the  grating  vector  orientation  and  the  angle  of  incidence.  They  can  be  characterized  as 
geometrical  conditions  since  they  involve  geometrical  requirements.  Condition  (4)  can  be 
characterized  as  amp/itucfe/phasecondition  since  it  involves  the  amplitudes  and  the  phases 
of  the  diffracted  plane  wave  components.  In  the  following  subsection  the  above  conditions 
are  discussed  and  analyzed  for  the  two  optic  axis  orientations  that  are  of  interest. 


A.  Geometrical  Conditions 

Al.  Optic  Axis  Perpendicular  to  the  Electrode  Surface 

The  case  of  the  optic  axis  perpendicular  to  the  electrode  surface  corresponds  to  the 
decoupled  case  where  both  TE  and  TM  modes  can  be  exist.  Since  the  optic  axis  is  along 
the  z  axis  there  is  a  cylindrical  symmetry  in  the  xy  plane.  Various  planar  slices  through 
the  wavevector  surfaces  are  shown  in  Fig.  4.  In  this  case  only  two  waves  (both  ordinary 
or  both  extraordinary)  exist  for  one  specific  mode.  This  b  the  reason  for  the  two  different 
propagation  constants  0rz>  shown  in  Figs.  4c  and  4d  for  a  TE  and  a 

TM  mode  respectively.  A  degenerate  case  of  both  a  TE  and  a  TM  mode  having  the  same 
propagation  constant  b  abo  possible  but  this  b  not  a  general  case  of  interest. 

It  b  straightforward  to  show,  using  the  wavevector  dbpersion  relations,  that  the 
normalbed  (with  respect  to  A:o)  V  components  of  the  diffracted  wavevectors  for  the  i-th 
diffracted  order  are 


pi) 


for  a  TE  incident  mode,  and 


±(/^«  - 

ngj 


1/3 
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(23) 

(24) 


where  the  +  sign  corresponds  to  backward-  and  —  sign  to  forward-diffracted  orders  and 
k,i  =  P'sinS  —  iK,  (in  this  case  due  to  the  unslanted  grating  K,  =  K  and  Ky  =  0). 
The  are  the  y  components  for  the  t-th  diffracted  order  of  a  p  polarization  (p  =  O,  E) 
incident  and  a  g  polarization  (9  =  O,  E)  diffracted.  Equations  (21),  (22)  and  (23),  (24) 
hold  for  both  0+,  O—  waves  and  E+,  E—  waves  respectively.  Consequently,  condition 
(1)  is  satisSed  since  both  0+,  O—  or  E+t  E—  waves  for  a  TE  or  a  TM  incident  mode 
respectively  lie  in  the  same  vertical  plane.  Using  Eqs.  (2l)-(24)  the  diffracted  angles  of  any 
diffracted  order  can  be  found  from  tanff**  =  These  angles  are  shown  in  Fig.  5 

for  the  0-order  and  the  t-order  forward-diffracted  waves.  In  the  case  of  a  TE  incident  mode 
(Fig.  5a)  both  plane  wave  components  in  the  directions  and  69^  are  ordinary  and 
have  the  same  propagation  constant  as  the  incident  mode.  In  the  directions  characterized 
by  the  angles  and  ,  there  are  two  extraordinary  waves  produced  by  the  OE  type 
diffraction.  These  two  extraordinary  waves  can  constitute  a  TM  mode  if  {ne /  /nof)P!rE 
is  a  propagation  constant  corresponding  to  a  TM  mode.  A  similar  situation  is  shown  in 
Fig.  5b  for  a  TM  mode.  The  two  ordinary  waves  along  the  direction  (produced  by 
the  EO  type  diffraction)  can  constitute  a  TE  mode  if  (no//ns/)^^^  is  a  propagation 
constant  that  corresponds  to  a  TE  mode. 

Due  to  the  cylindrical  symmetry  in  the  xy  plane  (shown  in  Fig.  4b),  it  is  straight¬ 
forward  to  prove  that  the  eigenvectors  and  the  eigenvalues  of  matrices  [Eq.  (Al), 
*  =  ft  <1  independent  of  the  angle  of  incidence  S.  The  i-th  diffracted  mode  has  the 
same  propagation  constant,  or  ,  for  TE  or  TM  mode  respectively  (Figs.  4c  and 
4d).  Thus,  Eqs.  (3)-(5)  are  satisSed  for  any  diffracted  order.  Consequently,  condition  (2) 
b  also  sat'isfied. 

The  last  of  the  geometrical  conditions  deab  with  the  Bragg  conditions.  The  Bragg 
condition  for  TE  or  TM  modes  can  be  derived  from  the  expressions  of  Ref.  58.  Both 
0+  and  O-  waves  satisfy  the  same  Bragg  condition.  In  addition,  E-\-  and  E-  share  the 
same  Bragg  condition.  For  a  TE-to-TE  or  a  TM-to-TM  mode  diffraction,  in  the  unslanted 
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electrode  configuration,  the  Bragg  condition  is  A,  =  Ao  «72;9' sin  ^ ,  where  q  =  TE,  TM, 
and  Ao  is  the  freespace  incident  wavelength.  The  geometrical  conditions  discussed  in  this 
section  are  extended  in  the  next  section  for  the  case  of  the  optic  axis  lying  in  the  plane  of 
the  electrodes. 


A2.  Optic  Axis  in  the  Plane  of  the  Electrode  Surface 

The  case  of  the  optic  axis  lying  in  the  electrode  plane  (xy  plane)  corresponds  to 
the  hybrid  mode  case  since  neither  TE  nor  TM  modes  can  exist  in  this  case.  The  plane 
intersections  of  the  wavevector  surfaces  corresponding  to  that  orientation  of  the  optic  axis 
are  shown  in  Fig.  6. 


Similarly  to  the  previously  case,  it  is  straightforward  to  show  that  the  normalized  y 
components  of  the  diffracted  wavevectors  are 

=*(^'’-*2*)'".  (25) 

I  B  B  _  ^ 


sin*  6f‘  - 

sin  26®*  + 

f 

=  i(«o/  “ 

”o/ 

c„  sin*  — 

e,„sin25**  -1- 

£„  cos’ 

1 

r 

where  again  the  +  sign  corresponds  to  backward-  and  the  -  sign  to  forward-diffracted 
wavevectors,  and  all  the  relative  permittivities  are  expressed  in  the  (z,  y,  z)  coordinate 
axis  system.  The  diffracted  angles  for  any  diffracted  order  can  be  found  from  the  equation 
tantff*  =  ±k,ifl^*.  These  angles  are  shown  in  Fig.  7.  In  the  EE  and  EO  type  of 
diffractions  the  computation  of  the  and  requires  the  solution  of  a  trigonometric 
equation. 
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In  the  general  case  there  are  four  distinct  propagation  directions  for  one  diffracted 
order  due  to  the  various  types  of  diffraction.  In  Fig.  7a  the  different  diffracted  angles  are 
shown  for  the  0-order  and  t-order  forward-diffracted  waves.  The  eight  diffracted  waves  are 
distributed  as  follows:  two  ordinary  waves  O-f-,  O—  along  the  direction  of  B9°  produced 
by  the  OO  type  diffraction,  two  extraordinary  waves,  J5+,  E—  along  the  direction  5** 
produced  by  the  EE  type  of  diff'raction,  two  ordinary  waves  0+,  0~  along  the  direction 
produced  by  the  EO  type  diffraction,  and  two  extraordinary  waves,  E+^  B-  along 
the  direction  *  produced  by  the  OE  type  of  diffraction.  However,  a  situation  like  the 
one  depicted  in  Fig.  7a  does  not  satisfy  geometrical  condition  (1)  since  for  a  mode  all  four 
waves,  0+,  0-,  E-\-,  and  E—  do  not  lie  in  the  same  vertical  plane.  This  problem  does 
not  appear  for  the  0-order  forward-diffracted  mode  since  as  is  shown  in  Fig.  7a  all  four 
waves  lie  in  the  same  vertical  plane.  There  is,  however,  a  solution  to  that  problem  if  the 
orientation  of  the  optic  axis  is  restricted  to  be  along  the  z  or  the  y  coordinate  axis.  In  this 
case  it  can  be  shown  that 

(M) 

In  addition  to  the  restriction  in  the  optic  axis  orientation  the  Bragg  condition  should  abo 
be  satisfied  in  order  for  Eq.  (29)  to  hold.  This  situation  is  shown  in  Fig.  8  for  the  two  optic 
axis  orientations.  It  is  important  to  note  the  symmetry  of  these  two  special  optic  axis 
orientations.  The  optic  axis  orientation  is  perpendicular  or  parallel  to  the  grating  vector. 

The  second  geometrical  condition  is  also  needed.  The  normalized  propagation  constant 
along  the  0-order  and  the  Bragg-order  forward-diffraction  directions  is  the  same  with  the 
normalized  constant  of  the  incident  mode,  The  condition  expressed  by  Eq.  (3)  should 
also  hold  to  guarantee  the  guidance  of  the  mode.  This  condition  for  the  incident  mode  can 
be  written  as 

= 0-  (30) 

It  is  straightforward,  but  tedious,  to  show  that  the  Bragg  diffracted  order  at  an  angle 
—S,  as  it  is  shown  in  Fig.  8,  will  satisfy  the  condition  (30)  for  ~S  and  fi'.  Consequently, 
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the  Bragg-order  forward-diffracted  waves  at  the  direction  =  —6  satisfy  the 

second  geometrical  condition.  The  third  geometrical  condition  is  concerned  with  the  Bragg 
conditions.  If  the  optic  axis  is  along  the  x  or  the  y  direction  and  the  Bragg  condition  is 
satisfied  for  the  ordinary  and  the  extraordinary  waves,  then  all  four  plane  wave  components 
of  the  incident  mode  will  satisfy  the  same  Bragg  condition  that  was  given  for  the  decoupled 
cases  where  This  can  be  seen  from  Fig.  8  and  can  be  found  analytically  using  the 

expressions  for  the  Bragg  conditions. 

The  rest  of  the  plane  waves  at  the  directions  of  the  angles  and  cannot 
constitute  a  mode  since  there  are  either  only  ordinary  or  only  extraordinary  components 
in  their  respective  directions,  and  the  modes  for  these  directions  must  be  hybrid  modes.  In 
addition,  all  the  diffracted  orders  except  the  0-order  and  the  Bragg-order  forward-diffracted 
modes  cannot  constitute  a  mode  since  they  do  not  satisfy  all  the  geometrical  conditions. 
The  modes  that  are  diffracted  in  these  directions  will  radiate  into  the  substrate  and  the 
superstrate,  and  consequently,  the  guided  mode  will  experience  radiation  losses  during  the 
diffraction  by  the  interdigitated-electrode  induced  phase  grating. 


B.  Amplitude  and  Phase  Matching  Conditions 

The  diffracted  plane  wave  components  of  the  four  incident  homogeneous  plane  waves 
(that  constitute  the  incident  mode)  can  be  written  in  the  form 

where  p  =  O,  E,  and  i  corresponds  to  the  »-order  forward-  or  backward-diffracted  wave. 
The  above  equation  can  be  derived  using  the  expressions  of  Ref.  58,  where  ordinary  and 
extraordinary  components  have  been  distinguished  and  the  exp(— factors  (s  is 
the  grating  thickness)  have  been  included  in  the  complex  T^ip±  (w  =  x,  y,  z)  wave 
amplitude  components.  Equation  (31)  includes  only  the  plane  wave  components  that  lie 
in  the  vertical  plane  that  make  an  angle  with  the  y  axis,  since  as  was  justified 
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previously,  all  four  plane  wave  components  0+,  O—,  £’+,  and  will  lie  in  that  plane 
if  the  geometrical  conditions  are  satisfied.  The  total  field  at  the  direction  of  5"  \b  Ei  = 
E^'*’  +  E°~  +Ef*  •+Ef~  .  These  equations  are  expressed  in  the  (z,  y,  t)  coordinate  axis 
system.  Now  it  is  convenient  to  express  the  field  components  in  the  waveguide  coordinate 
axis  system  (this  system  b  shofwn  in  Figs.  5  and  7).  It  b  worth  mentioning  that  the 
waveguide  coordinate  system  b  defined  as  a  function  of  the  diffracted  angle.  For  the  t- 
order  forward-diffracted  wave  the  transformation  between  the  waveguide  coordinate  system 
{^wiy  y»»»  *••)  coordinate  axis  system  (i,  y,  z)  b  a  rotation  by  and  the 

translation  along  the  z  axis  has  been  suppressed  because  does  it  not  affect  the  results. 
Using  that  transformation  the  the  diffracted  fields  are: 

where  (u  =:  z,  y,  z)  are  the  complex  amplitude  components  expressed  in  the 

(2«.«  Vwi,  z«<)  coordinate  axb  system. 

Bl.  Optic  Axis  Perpendicular  to  the  Electrode  Surface 

Using  the  properties  of  the  ordinary  and  the  extraordinary  waves  and  the  geometrical 
conditions  it  b  straightforward  to  show  for  TE  modes  that 

=  y,iT,^io±  exp[-y^o(±*.,,o*-<  +  (33) 

and  for  the  TM  modes  that 

Elf  =  -f-7;^.E*i.iexpl-yA:o(±Ar,.,£Z.,  -l-^'z.,)].  (34) 

Consequently,  the  »-order  diffracted  waves  of  a  TE  incident  mode  will  constitute  a  TE 
mode  and  i-order  diffracted  waves  of  a  TM  incident  mode  will  constitute  a  TM  mode.  It  b 
convenient  to  compare  the  diffracted  electric  fields  with  the  incident  electric  fields,  which 
arc  given  by  Eq.  (A5)  for  a  TE  mode  and  by  Eqs.  (A9)  and  (AlO)  for  a  TM  incident  mode. 
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Thus,  the  complex  amplitude  ratio  coefficients  are  defined  for  the  TE  modes  (ordinary 
waves)  as 


TywiO  + 


TywiO- 


fio-i-  =  l/.o+lexp(-;Vio  +  )  = 

/<o-  =  l/io-|exp(-jV.o-)  =  p  , 

ro- 

and  for  the  TM  modes  (extraordinary  waves) 


/\  f  I _ /  -r  \  + 

iS+  ~|/»£+l®^P(  J^0^»£+)  tfiilg-  \F  ~  (L  Ini  \  r 

At- =I/<£-I«xp'-Ap,j.)= 


(35) 

(36) 


(37) 

(38) 


I^^oi)Fb- 

The  /,o+  ,  /,o-  I  /.£+  »  /.£-  coefficients  express  ratios  between  incident  and  diffracted 

complex  amplitude  components.  Using  the  above  defined  coefficients  the  t-order  diffracted 
total  electric  field  can  be  written  for  the  TE  modes  as 


FyM,xj  —  [/iO+  ^o  +  <*o+  (*•»■  )  +  fio~Fo-ao.  (x,.)]a, (£.,), 


(39) 


where  ao+  (x,<)  and  Cq-  (x.<)  were  defined  previously  (as  a  function  of  z.)  and  o,  (z,,)  = 
exp(— Similar  expressions  to  those  of  Eqs.  (A6)  and  (A7)  can  be  written  for 
the  total  magnetic  field  that  corresponds  to  the  i-order.  Now  for  the  TM  modes  the  total 
electric  field  components  for  the  t-th  order  are 


Famif  — 


P' 


^mmwj 


(/<£+ ^£ +  «£+ (X.<)  +  As  - -  0£ -(*•<  )l“i  (*•<). 


E.nH  =  -^^^[A£  +  f'£4«£+(x.<)- 

^O! 


(40) 

(41) 


The  complex  coefficients  fio^  •  Uo  - » /•£+ » /•£-  influence  the  phase  and  the  amplitude 
relationship  of  the  plane  wave  components  of  the  diffracted  mode.  The  mode,  however, 
should  satisfy  the  characteristic  mode  amplitude  and  phase  relationship  for  its  plane  wave 
components.  Thus,  only  the  in-phase  waves  with  the  correct  amplitude  ratios  as  given  by 
Eq.  (39)  for  the  TE  modes  and  by  Eqs.  (40)  and  (41)  for  the  TM  modes  will  constitute 
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the  i-th  diffracted  mode.  These  phase  and  amplitude  parts  of  the  electric  fields  for  the  TE 


modes  axe 


^ywif  Ao  +  PT  fi  [y*  J  4  +  flo  4  t )  d"  i  )  j 


and  for  the  TM  modes  are 


^•wi/  —  [As+  ■f£  +  <*s+  (*•<)  +  As-  Fb-Ce-  (x*<)]oi(z,<),  (43) 

*»«»  / 

E.,>,  =  -  (i..)  -  r.-  (r.,)K(».,),  (44) 

"o/ 

where  p,  =  min{l,  max{(lAr-  |/|Ab4  |)cos(^iB+  -  (Pis-),  0}},  and  g  =  O,  E  {O  corre¬ 
sponds  to  TE  modes  and  E  corresponds  to  TM  modes).  The  coefficients  po  and  pg  satisfy 
the  inequality  0  <  po ,  Pa  <1,  and  they  represent  the  radiation  losses  of  the  mode  due  to 
the  phase  and  amplitude  mismatches.  For  this  decoupled  optic  axis  orientation,  however, 
due  to  the  symmetry  of  the  problem  Ao+  —  Ao-  for  the  TE  modes  and  As+  As-  for 
the  TM  modes.  Consequently,  po  1  and  Pg  -  I  which  means  that  in  this  decoupled 
case  there  is  negligible  radiation  loss  due  to  the  amplitude  and  the  phase  mismatches. 


B2.  Optic  Axis  in  the  Plane  of  the  Electrode  Surface 


In  this  case  the  modes  are  hybrid  and  the  general  equation  (32)  must  be  used.  Again 
coefficients  Ao-f  >  Ao  - 1  Aa-t-  >  and  As-  can  be  defined  as  follows 


A>±  lA>±lc^rp(  JPipi)  ~ 


T 

•‘■mwipfpm 
0 _ 


-f'p± 


- p  a7  „  *  v’-; 

Za 

*— /r 

where  p-^  O,  E.  Using  a  similar  procedure  to  the  one  above  for  the  TE  and  TM  modes, 
the  following  electric  field  components,  that  satisfy  the  phase  and  amplitude  conditions, 
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can  be  written 


^MWif 

Eywif 


^  *-**^*^'- (x..)  -  F«-a,,(x,,)la.(e.O. 

=  o.e 

/<o+p[^’,+  a«+(a:,.)  +  F,.a,_(x,J]a<(e,.), 
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+  F,_a,_  (x«<)]ai(«-<). 
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where 


(46) 

(47) 

(48) 


/>  =  min{l,  mai{}^^cos(^iO-  -(Pio*),  0},  mai{}4s^  cos(ic>,n.  -Viot),  0), 
l/tO+l  l/»o  + 


-V><o  +  ),  0)}. 
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(49) 


Again  the  coeflicient  p  is  related  to  the  radiation  losses  due  to  the  phase  and  amplitude 
mismatches  of  the  plane  waves  and  satisfies  the  inequality  0  <  p  <  1.  Similar  expressions 
to  Eqs,  (46)-(48)  can  be  written  for  the  magnetic  fields.  Figure  9  demonstrates  the  general 
relation  of  the  complex  coefficients  /<©+ ,  /<o-  »/<£+»  a^nd  /i£_ ,  in  the  complex  plane. 


If  p  =  0,  which  is  a  possible  situation,  then  the  diffracted  mode  will  be  coupled  to 
radiation  modes  only.  In  the  following  section  some  comments  on  the  diffracted  mode 
efficiencies  are  discussed. 


VII.  Mode  Diffraction  EfficiencicB 

In  order  to  determine  the  efficiencies  of  the  diffracted  modes  it  is  important  to  compute 
the  power  per  unit  area,  along  the  x^i  and  coordinates,  that  propagates  in  the  x., 
direction.  The  power  per  unit  area  can  be  found  at  any  point  using  the  components  of  the 
Poynting  vector.  Then,  the  power  per  unit  area  that  b  propagating  in  the  z^i  direction 
can  be  found  by  a  simple  integration 

PawiS  ~  f  ^»»t/ (*»*)^^»* »  (^®) 

Jo 
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where  P,^ij  is  the  power  per  unit  area  in  the  direction  and  5,,i/  is  the  component 
of  the  Poynting  vector  as  given  in  Appendix  C.  Both  P,«i/  and  correspond  to  the 
film  region  only.  The  power  in  the  cover  and  the  substrate  can  be  computed  similarly,  but 
since  it  was  assumed  that  the  evanescent  fields  in  the  cover  and  the  substrate  regions  follow 
the  same  diffraction  characteristics  as  the  corresponding  plane  waves,  the  expressions  for 
the  film  region  power  will  be  generalized  for  the  total  power  propagating  along  the 
direction  of  the  waveguide. 


The  power  per  unit  area  along  the  direction  can  be  found  using  Eq.  (50).  The 

efficiency  of  a  diffracted  mode,  DEi ,  can  be  defined  as 

,  cos  £*• 


DEi  = 


P»»*7  _  a 


(51) 


n  r  +  I  r  * 

Pmw/ 

where  P,,/  is  the  incident  power  per  unit  area  and  q  =  O,  E.  If  q  =  O,  Eq.  (51)  gives  the 
diffraction  efficiency  of  a  pure  TE  mode  while  if  q  =  E  it  gives  the  diffraction  efficiency  of 
a  TM  mode.  For  the  hybrid  mode  case  69^  =  =  Si  and  q  =  O. 


Vin.  Diffraction  bv  Slanted  Electrode  Configurations 

In  this  section  the  diffraction  of  guided  waves  in  anisotropic  dielectric  slab  waveguides 
by  slanted-interdigitated-electrode  induced  phase  gratings  will  be  discussed.  Again  the 
mode  propagation  conditions  described  in  section  VI  should  hold.  The  azaplitude/pbase 
conditions  do  not  change  in  the  case  of  slanted  electrodes  and  the  analysis  presented  in 
subsection  VI.B  is  still  valid  under  the  same  given  conditions  of  section  I.  There  are, 
however,  some  differences  in  the  geometrical  conditions  that  will  be  identified  for  the  two 
optic  axis  orientations  of  interest. 


A.  Optic  Axis  Perpendicular  to  the  Electrode  Surface 

Due  to  the  cylindrical  symmetry  of  the  wavevector  surfaces  (along  the  z  axis),  all 
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the  geometrical  conditions  are  still  satisfied,  and  Eqs.  (21)-(24)  are  valid.  In  this  case, 
however.  A",  =  Asin^,  where  ^  is  the  slant  angle  of  the  induced  grating.  Again,  two 
distinct  decoupled  mode  types  can  be  identified,  TE  and  TM  modes.  In  Fig.  10  the  zy 
plane  wavevector  intersections  are  shown  for  a  TE  mode  (Fig.  10a)  and  for  a  TM  mode 
(Fig.  10b).  The  cylindrical  symmetry  for  both  modes  types  does  not  affect  the  geometrical 
conditions  of  the  problem.  The  Bragg  conditions  are  given  by  A,  =  Aot72/?'  cos(^  -  d), 
where  q  =  O,  E,  for  the  two  ordinary  or  the  two  extraordinarj'  waves  of  the  incident  mode. 


B.  Optic  Axis  in  the  Plane  of  the  Electrode  Surface 

This  is  the  case  of  a  hybrid  incident  mode  where  both  ordinary  and  extraordinary 
homogeneous  plane  waves  are  incident.  Multiple  Bragg  conditions  and  difi'raction  direc¬ 
tions  are  possible  similar  to  those  shown  in  Fig.  7.  Equations  (25)'-(28)  hold  again  if 
Ac  =  Asin^.  In  order  to  satisfy  the  geometrical  conditions,  a  restriction  in  the  orienta¬ 
tion  of  the  optic  axis  is  required.  It  can  be  easily  shown  that  when  the  optic  axis  is  parallel 
or  perpendicular  to  the  grating  vector  orientation  and  the  Bragg  condition  b  satbfied  then 

k^o  ^  ^  ^^(2^  -  6 -b),  (52) 

and  the  four  homogeneous  plane  wave  components  lie  in  the  same  vertical  plane  at  an 
angle  69°  =69^  =  2^  —  d  —  x  and  the  same  propagation  constant  fi'.  The  sign  in 
Eq.  (52)  should  be  chosen  such  that  k°P  =  <  0  (propagation  along  the  -y  direction). 

Consequently,  when  the  optic  axb  b  parallel  or  perpendicular  to  the  grating  vector  and  at 
Bragg  incidence  the  first  geometrical  condition  b  satbfied.  The  second  condition  b  abo 
satbfied  since  Eq.  (30)  holds  for  2^-6— r  abo,  as  can  be  shown  using  the  same  propagation 
constant  0’.  Finally,  the  third  geometrical  condition  b  abo  satisfied  and  all  the  ordinary 
and  extraordinary  incident  plane  waves  share  the  same  Bragg  condition  given  as  before 
(for  =  /(?').  Plane  intersections  of  the  wavevector  surfaces  are  shown  in  Fig.  11  for  the 
optic  axis  parallel  to  the  grating  vector  (Fig.  11a)  and  for  the  optic  axis  perpendicular  to 
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the  grating  vector  (Fig.  11b). 


DC.  Multiple  Slanted  Electrode  Configuratione 

It  is  relatively  straightforward  to  generalize  the  above  results  for  the  case  of  multi* 
pie  cascaded  gratings.  The  ampUtude/phiLse  conditions  remain  unaffected  if  the  multiple 
cascaded  anisotropic  grating  diffraction  analysis  described  in  Ref.  58  is  used.  The  geomet¬ 
rical  conditions  in  the  case  of  the  optic  axis  being  perpendicular  to  the  electrode  surface 
are  again  satisfied  due  to  the  cylindrical  symmetry  of  the  wavevector  surfaces  (along  the 
t  axis).  Consequently,  any  cascaded  configuration  of  electrodes  can  be  treated  without 
restrictions  if  the  optic  axis  is  perpendicular  to  the  electrode  surface.  In  addition,  if  the 
optic  axis  lies  in  the  electrode  plane  then  the  grating  vector  should  be  parallel  or  perpen¬ 
dicular  to  the  electrode  surface  and  if  4>t  and  i  are  the  slant  angles  of  two  consequetive 
gratings  then  =  0,  or  x/2.  An  additional  constraint  is  that  the  multiple  type 

Bragg  conditions  must  hold  in  order  to  satisfy  all  the  geometrical  conditions.  Summarizing, 
the  hybrid  mode  case  is  much  more  restrictive  than  decoupled  cases  due  to  the  differing 
diffraction  characteristics  of  the  ordinary  and  extraordinary  waves. 


X.  Lithium  Niobatc  Example 

In  this  section  two  examples  of  guided  wave  anisotropic  diffraction  analyses  are  de¬ 
scribed.  Those  examples  are  based  on  two  applications  of  lithium  niobate  interdigitated 
electrode  devices.  These  applications  are:  a  Givens  rotation  device  and  a  herringbone  mul¬ 
tiplier.  In  both  of  these  «xamples  the  optic  axis  orientation  is  chosen  to  be  perpendicular 
to  the  electrode  surface.  This  configuration,  as  will  be  justified  later,  is  the  most  efficient 
if  the  waveguide  is  designed  to  support  only  the  TMo  mode.  This  can  be  achieved  if  the 
change  in  the  extraordinary  film  index  is  greater  than  the  change  of  the  ordinary  film  index 
comparing  these  indices  with  the  corresponding  substrate  indices.  In  all  the  examples  the 
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film  and  substrate  material  is  lithium  niobate  while  the  cover  is  silicon  dioxide  which  serves 
as  a  buffer  layer  between  the  electrodes  and  the  waveguide.  The  superstrate  is  air.  The 
operating  freespace  wavelength  is  taken  to  be  1.0/im  and  the  refractive  indices  (at  that 
wavelength)  arc  =  1.51,  no/  =  2.2512,  fisf  =  2.1894,  no,  =  2.2362,  and  ng,  =  2.1594. 
These  indices  correspond  to  a  titanium  diffusion  formed  waveguide.  The  film  thickness  is 
taken  to  be  l.O/zm.  Solving  the  waveguiding  problem  using  Eq.  (3)  for  T£  modes  or  TM 
modes  (due  to  the  decoupled  orientation  of  the  optic  axis  that  equation  may  be  greatly 
simlified  as  was  described  in  section  II.A).  For  the  chosen  film  thickness  only  TEo  and 
TMo  modes  can  propagate  with  normalized  propagation  constants  of  =  2.2367  and 
0^^  =  2.1671.  The  normalized  components  of  the  ordinary  or  extraordinary  wavevec- 
tors  of  the  plane  wave  mode  components  can  be  found  using  Eq.  (B6)  or  (B7)  respectively 
of  Appendix  B.  The  induced  relative  permittivity  tensor  under  the  electrodes  (inside  the 
waveguide)  depends  on  the  chosen  periodicity  of  the  grating  (or  equivalently  the  electrode 
spacing),  the  thickness  of  the  buffer  layer,  the  ratio  of  the  electrode  width  to  electrode 
spacing  (which  in  all  the  examples  is  tu/A,  =  0.50),  and  the  ratio  of  the  electrode  thick* 
ness  to  the  electrode  spacing.  In  addition,  the  electrostatic  relative  permittivities  are 
€e  ~  3.90  (SiOa),  and  e,,4  =  —  84.142  and  e.,4  s  28.122  (lithium  niobate).  The 

relative  permittivity  of  the  air  is  assumed  to  be  unity.  The  small  difference  in  the  rela* 
tive  permittivities  of  the  film  and  the  substrate  regions  is  negligible  in  the  analysis  of  the 
electrostatic  problem. 


A.  Givens  Rotation  Device  Example 

The  Givens  rotation  operation  occupies  a  central  role  in  linear  algebraic  signal  pro¬ 
cessing.  An  integrated-optical  coherent  implementation  of  an  elementary  rotation  matrix 
device,  can  be  constructed  based  on  a  thick  grating  induced  by  interdigitated-electrodes. 
Consequently,  the  same  device  that  has  been  described  in  the  previous  subsection  can  be 
used  as  a  Givens  rotation  device^’  to  implement  the  Givens  orthogonalization.  However, 
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in  this  case  the  Amplitude  of  the  transmitted  and  diffracted  fields  is  needed  for  correct  op¬ 
eration.  The  transmitted  field  must  vary  as  8in(if,V)  and  the  diffracted  field  as  cos(K,V) 
where  K,  is  a  constant  and  V  is  the  voltage  applied  to  the  in  ter  digitated  electrodes. 

The  Givens  rotation  device  is  shown  schematically  in  Fig.  12.  The  periodicity  of  the 
induced  grating  (and  consequently  the  electrode  spacing)  is  such  that  it  satisfies  the  first 
Bragg  condition  for  a  TM©  incident  mode.  Using  the  normalized  propagation  constant 
=  2.1671  for  an  angle  of  incidence  S  =  5®,  a  Bragg  periodicity  of  Aru  =  Ao  = 
2.6473/im  can  be  computed.  For  a  film  thickness  of  l.O/im  a  TEo  mode  can  also  propagate. 
However,  the  coupling  between  TMo  -  TEo  is  negligible  since  the  deviation  from  the  Bragg 
periodicity  can  easily  be  computed  from  the  equation, 

0?,  =  ^  (|  -  ')•  (“) 

and  is  found  to  be  0.25/im  which  is  very  large  for  “thick"  grating  diffraction.  The  buffer 
layer  thickness  is  d,  =  0.25pm,  and  the  electrode  thickness  is  2t  =  0.10pm.  Using  the 
analysis  described  in  Ref.  57  the  relative  permittivity  tensor  can  be  calculated  and  is 
eji  =  -0.39  X  10- ‘  V/pm,  =  -1.14  x  10" ‘  V/pm,  c,,  =  -0.15  x  lO"*  V/pm, 
Cii  =  0.69  X  10"*  V/pm,  and  the  ss  0.24  computed  with  averaging  along  the  z 
direction  of  the  electric  field  without  any  optical  weighting. 

The  0-order  and  the  -l-l-order  forward-diffracted  electric  fields  are  shown  in  Fig.  13 
as  a  function  of  the  applied  voltage  normalized  with  respect  to  the  incident  electric  field. 
The  dotted  lines  correspond  to  rigorous  calculations  using  the  analysis  described  in  this 
chapter,  while  the  continuous  lines  correspond  to  a  cosine,  cos(if,  V),  (Fig.  13a)  or  a  sine, 
8in(if,V),  (Fig.  13b)  best  fit  on  the  rigorously  computed  results.  The  best  fit  value  for 
K,  was  0.123  rad/volt.  The  agreement  between  the  rigorous  curves  (dotted  lines)  and  the 
cosine  and  sine  approximations  (continuous  curves)  is  excellent.  This  good  agreement  is 
necessary  for  the  correct  device  operation. 
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B.  Herringbone  Multiplier  Example 

Similarly  to  the  previously  described  case,  a  Z-cut  (optic  axis  along  the  z  direction) 
slab  is  used  for  the  design  of  the  herringbone  multiplier.  It  is  taken  that  the  TMo  mode  is 
propagating.  The  angle  of  incidence  is  5  =  —5®  (Fig.  la).  The  TMo  mode  consists  of  two 
extraordinary  polarized  plane  waves  with  parameters  that  can  be  computed  by  £q.  (B7), 
where  again  =  2.1671.  The  slant  angle  of  the  first  grating,  <t>i »  can  be  computed 

using  Ref.  58  for  t'l  —  =  1  to  be  =  87.5®.  Another  solution  of  the  same  equation 

is  =  —2.5®  which  corresponds  to  a  backward  diffracting  grating  and  consequently  is 
not  practical  for  the  design  of  the  multiplier.  The  slant  angle  of  the  second  grating  is 
^3  =  92.5®  (again  using  the  equations  of  Ref.  58).  The  common  periodicity  A  of  both 
gratings  c2Ln  be  calculated  to  be  5.Znm.  The  thicknesses  of  both  gratings  were  chosen 
to  be  equal,  Si  =  S3  =  1000/xm.  The  top-view  of  the  multiplier  and  the  corresponding 
wavevector  diagram  are  shown  in  Figs.  14a  and  14b  respectively.  For  the  electrostatic 
analysis  the  ratio  of  the  electrode  gap  to  the  electrode  spacing  is  t/Kt  =  0.5  (where  A. 
is  the  electrode  spacing).  The  buffer  layer  is  assumed  to  be  O.lSfxm  and  the  electrode 
thickness  0.25pm.  Using  the  electrostatic  analysis  of  Ref.  57  the  fundamental  component 
of  the  electric  potential  Qi  is  found  to  be  0.3146  volts.  The  induced  permittivities  are 
e..k  =  (ffii*  cos[Ke  •  f)  +  e334  8in(J?  •  0](V*/A«fc), 

=  l-Ciifc  cos(Ke  •  f)  +  e,3k  Bin(K-f)](V’fc/A.k), 
e..k  =ffs9fc{6in(^i  •OK^fc/A.k),  ^  ^ 

,  (54) 

=  Cx3fclcos(K<  •  OK'^kM.k), 

e.,k  =  e,sklcos(K<  •f)l(V'fc/A.k), 
e,,k  =  C39»lcos(Kt  •f)](V'fc/A.fc), 

where  A;  =  1,  2,  ffm  =  0.32  x  10~*  V/fitn,  Cm  —  Caai  —  faaa  —  ~0.19  x 

10* *  V/pm,  Cssi  =  ejsa  =  -0.54  x  10"’  V/pm,  Cxai  =  ^laa  =  -0.72  x  10"*  F/pm, 
ffisi  =  Cl sa  =  0.32x10"*  V/pm,  Caai  = -Sjsa  = -0.14  x  10"*  V/pm,  pj*  =0.43,  and  Fj, 
F3  are  the  applied  voltages.  In  order  to  operate  the  multiplier  it  is  a  requirement  to  have 


an  output  which  is  proportional  to  the  product  of  the  two  applied  voltages  V^V^.  Thus, 
the  desired  diffraction  efficiency  is  of  the  form  DE  jss  if,  Vi  V3 . 

The  performance  of  the  multiplier  has  been  evaluated  with  a  numerical  example. 
In  order  to  achieve  better  performance  a  set  of  bias  voltages  have  been  added.  The  bias 
voltages  V^i ,  Via  ^ere  chosen  10.8  volts  in  order  to  operate  as  close  as  possible  to  the  linear 
response  of  the  diffraction  curves  of  each  grating.  The  input  voltages  V, ,  Vj  were  varied 
around  the  bias  voltages  in  the  range  from  8.8  volts  to  12.8  volts.  The  -2-order  forward- 
diffracted  mode  efficiency  is  shown  in  Fig.  15  as  a  function  of  the  voltage  V  =  Vi  =  V3, 
along  with  the  predicted  efficiency  from  the  equation  DEpr  =  K,  Vi  V3  (where  K,  has  been 
computed  around  the  bias  point  to  be  (0.25/10.8^)  l/volt’).  The  error  as  a  function  of  the 
applied  input  voltage  is  also  shown  in  Fig.  15.  The  percentage  error  has  been  calculated 
using  the  formula  PE  =  {DEp^  —  DE)  x  100 fDE  where  DE  is  the  rigorously  computed 
efficiency  and  DEpr  is  the  predicted  efficiency  using  the  equation  DEpr  =  if,  V1V3.  It  is 
observed  that  the  error  is  in  the  order  of  ±8%  for  the  chosen  range  of  variation  of  the 
input  voltages  around  the  bias  points.  One  potential  problem  of  the  proposed  multiplier 
is  the  limited  dynamic  range  if  larger  voltages  are  necessary.  Consequently,  compensating 
electronics  may  be  necessary. 


XI.  Hybrid  Mode  Example 

In  this  hybrid  mode  case  the  lithium  niobate  crystal  is  assumed  to  be  Y-cut  and  X 
propagating.  Consequently,  the  optic  axis  u  along  the  x  axis  of  the  coordinate  axis  system 
that  is  used  in  this  analysis.  The  refractive  indices  of  the  cover,  film  and  substrate  are  n,  = 
1.51,  nos  -  2.254,  n*/  =  2.174,  no,  =  2.234,  and  n*.  =  2.154,  at  an  operatmg  freespace 
wavelength  of  Ao  =  l.Ofi/im.  The  thickness  of  the  film  region  b  again  assumed  to  be  l.O/xm. 
According  to  Eq.  (6),  in  order  to  have  four  homogeneous  plane  wave  components  the  angle 
of  incidence  6  should  be  greater  than  a  critical  angle  which  for  the  given  parameters  of 
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the  waveguide  b  60.7” .  In  this  analysis  the  angle  of  incidence  is  chosen  to  be  f  =  62”  and 
therefore  is  greater  than  the  critical  angle.  Using  the  above  defined  parameters  and  the 
analysis  of  Refs.  59  and  60  (or  Eq.  (8)|,  two  acceptable  values  of  the  normalized  propagation 
constant  0'  =  2.217,  and  /?'  =  2.235  are  calculated.  From  these  values  only  /?'  =:  2.235 
corresponds  to  four  homogeneous  plane  wave  mode  components.  Using  Eqs.  (B6)  and  (B7) 
of  Appendix  B  all  the  plane  wave  parameters  can  be  computed.  From  the  Bragg  condition 
A  =  0.27/im.  This  periodicity  is  very  small  for  practical  applications  due  to  the  large 
angle  of  incidence  but  the  results  of  this  example  analysis  are  only  for  demonstration  of 
the  hybrid  mode  effects.  Using  a  buffer  layer  of  O.OSAim  and  f/A,  =  0.5,  a  fundamental 
electric  potential  harmonic  of  Qi  =  0.054  volts  was  computed.  The  induced  permittivity 
tensor  is 

e..  =  en[cos(ifi))(V/A.), 

=  l-Caa  cos(Kz)  -  C95  sin{Ki)](V/A«), 
e,,  s=  (caa  cos(Kx)  -  ess  sin(/fi))(V/A.),  (55) 

e„  «  e„  =  0, 

=  €,slsin(Kz)l(V/A,), 

where  Cn  =  0.61  x  10“  *  V/fim,  Caa  =  0.22  x  10"^  V/nm,  ess  =  0.85  x  10"*  V/fim, 
eis  =  —0.38  X  10“  *  V/nm,  and  pj*  or  0.5.  The  grating  thickness  was  assiimed  to  be 
2000(im,  and  the  applied  voltage  was  varied  between  0  and  50  volts.  The  diffraction 
efficiencies  of  the  0+,  0-,  E+,  and  E-  waves  are  shown  in  Fig.  16  and  their  corresponding 
phases  in  Fig.  17  (for  the  border  and  the  Bragg-order)  as  a  function  of  the  applied  on  the 
electrodes  voltage.  The  resulting  mode  efficiencies  are  shown  in  Fig.  18  for  the  0-order 
(dashed  line)  and  the  Bragg-order  (continuous  line).  It  is  shown  that  the  diffracted  mode 
efficiencies  in  this  case  appear  to  have  significant  radiation  losses  due  to  the  phase  and 
amplitude  mismatch  between  the  ordinary  and  the  extraordinary  waves. 
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XII.  Comparison  with  Experimental  Results 

In  this  section  the  results  of  the  analysis  presented  are  compared  with  the  experimental 
results  of  Ref.  18.  The  geometrical  configuration  is  that  of  Fig.  1.  The  electrodes  are 
unslanted  with  respect  to  the  (z,  y,  z)  coordinate  system  and  the  optic  axis  orientation 
is  along  the  z  axis.  Before  starting  the  analysis  of  that  configuration  it  is  worthwhile 
to  comment  on  this  optic  axis  orientation.  According  to  the  results  of  Refs.  59  and  60 
for  this  optic  axis  orientation,  the  propagating  modes  at  normal  incidence  to  the  grating, 
are  characterized  as  critically  stable.  Slight  divergence  in  the  angle  of  incidence  from 
propagation  along  the  — y  axis  changes  the  mode  from  pure  TEo  to  a  leaky  mode.  In  this 
leaky  mode  region,  due  to  the  negative  crystal  birefringence,  the  ordinary  component  of 
the  hybrid  mode  leaks  into  the  substrate. 

For  the  diffraction  process  it  is  necessary  to  use  angles  of  incidence  away  from  the 
normal.  For  this  propagation  direction  the  modes  are  leaky  guided  modes  and  power  will 
be  lost  as  it  travels  down  the  waveguide  even  if  the  grating  is  not  present.  Additional  losses 
are  also  present  due  to  the  diffraction  of  the  hybrid  mode  as  described  previously.  In  order 
to  study  the  above  experimental  configuration  the  hybrid  character  of  the  incident  mode 
will  be  suppressed  due  to  the  very  small  angle  of  incidence  (5  =  0.62”).  Consequently,  the 
leaky  character  of  the  mode  will  also  be  neglected. 

For  the  electrostatic  analysis  the  ratio  tfk,  is  0.5,  the  electrode  thickness  is  0.5/xm, 
and  the  buffer  layer  (SiOa)  thickness  is  0.15/xm.  The  grating  periodicity  A  is  13.33  ^m. 
The  low-frequency  relative  permittivities  of  all  the  materials  are  the  same  as  in  section  5.10. 
The  fundamental  harmonic  of  the  electric  potential  isQi  =  0.54  volts.  Using  that  value  of 
Qi  a  fundamental  electric  field  harmonic  of  Ei  —  0A22{Vg  /I)  is  computed,  where  Vg  b 
the  voltage  difference  between  two  consequetive  electrodes.  If  Engan’s  analysis  is  used’*  ** 
then  the  corresponding  fundamental  harmonic  isEi  ^  0.847(V](  /^).  As  a  result,  neglecting 
the  electrode  thickness  and  the  buffer  layer  introduces  an  error  in  the  electrostatic  field 
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computation  of  about  50%.  Similarly,  the  error  in  the  capacitance  computation  Lb  of  the 
same  order.  The  operating  freespace  wavelength  is  0.633/im  and  the  refractive  indices  of  the 
film  are  taken  to  be  no/  =  2.2885  and  n^/  =  2.2014.  At  normal  incidence  on  the  grating 
the  incident  mode  would  be  a  TEo  mode.  For  the  given  optic  axis  orientation  the  TEq  mode 
corresponds  to  two  extraordinary  incident  waves.  The  normalized  propagation  constant  is 
found  using  the  Bragg  condition  given  in  Ref.  18  and  is  2.194.  The  parameters  of  the  two 
extraordinary  waves  can  be  computed  from  Eq.  (11.7).  The  induced  relative  permittivity 
tensor  is  given  by  Eqs.  (55)  where  Cn  =  0.69x10”*  V//im,  C33  =  0.25x10“*  V//xm,  Ssj  = 
0.97  X  10“  *  V//im,  Sis  =  —0.43  x  10“  *  V//im,  and  g]*  =■  0.51.  The  Bragg  diffracted  mode 
efficiency  along  with  the  experimental  results  are  shown  in  Fig.  19.  An  important  difference 
between  the  experimental  and  the  theoretical  results  is  that  experimental  efficiencies  are 
normalized  with  respect  to  the  output  power  while  the  theoretical  results  are  normalized 
with  respect  to  the  input  power.  This  is  the  reason  that  the  losses  do  not  appear  in  the 
experimental  data  even  if  it  is  known  that  the  incident  mode  b  slightly  leaky. 


Xin.  Discussion  and  Summary 

It  has  been  shown  that  all  optic  axis  orientations  supporting  only  hybrid  modes  in¬ 
troduce  radiation  losses  after  the  diffraction  process  due  to  the  geometrical  and  ampli¬ 
tude/phase  mismatches  of  the  ordinary  and  extraordinary  waves.  In  some  cases  (leaky 
guided  modes)  thb  b  true  for  the  mode  even  without  the  diffraction  by  the  interdigitated- 
electrode  induced  grating.**'  **  In  thb  paper  the  most  efficient  optic  axb  orientation  b 
identified  to  be  the  orientation  corresponding  to  the  x  axb  direction  (perpendicular  to  the 
electrode  surface).  This  decoupled  optic  axb  orientation  guarantees  negligible  radiation 
losses  due  to  the  different  mismatch  types.  For  thb  preferred  orientation  of  the  optic  axb 
the  propagation  of  only  the  TMo  mode  can  be  achieved  (for  negative  birefringent  lithium 
niobate)  since  its  cutoff  thickness  b  smaller  than  that  for  TEo  mode.  Consequently,  the 
large  electrooptic  coefficient  can  also  be  utilized.  An  additional  advantage  b  the  increase 
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of  the  z  component  of  the  electrostatic  field  (that  affects  the  index  difference  due  to  the 
largest  electrooptic  coefficient  in  lithium  niobate).  As  a  result  the  orientation  of  the  optic 
axis  perpendicular  to  the  electrode  surface  should  be  preferable. 

In  the  case  of  the  hybrid  modes  the  resulting  diffraction  efficiencies  are  affected  by  radi¬ 
ation  losses.  This  is  due  to  the  various  mismatches  of  the  O-f,  O— ,  £-f-,  and  E—  waves  and 
especially  the  mismatch  between  the  ordinary  and  the  extraordinary  components.  How¬ 
ever,  the  analysis  neglects  the  coupling  between  the  different  guided  or  radiation  modes. 
This  coupling,  in  some  cases,  can  act  cooperatively  with  the  diffraction  process  resulting 
in  decreased  radiation  losses.  Consequently,  the  above  hybrid  mode  results  can  be  consid¬ 
ered  as  a  guideline  in  the  design  of  hybrid  mode  interdigitated-electrode  devices.  However, 
propagation  along  directions  where  leaky  modes  are  permittable  should  be  avoided. 

In  summary,  an  analysis  for  treating  the  diffraction  of  guided-waves  in  uniaxial  an¬ 
isotropic  slab  waveguides  has  been  developed  for  the  first  time.  The  analysis  is  based  on 
the  three-dimensional  anisotropic  rigorous  coupled-wave  diffraction  analysis  and  on  the 
decomposition  of  a  guided  mode  in  a  uniaxial  anisotropic  slab  waveguide  into  foiir  ho¬ 
mogeneous  plane  wave  components.  The  latter  decomposition  restricts  the  types  of  the 
guided  modes  that  can  be  treated  to  pure  guided  modes  which  are  decomposable  into  four 
homogeneous  plane  wave  components.  However,  the  analysis  for  the  first  time  takes  into 
account  the  anisotropic  properties  of  all  media,  an  arbitrary  number  of  diffracted  orders, 
and  the  noncollinear  character  of  the  diffraction  (where  the  incident  mode  and  the  grating 
vectors  are  not  collinear).  Example  analyses  have  been  presented.  These  include  a  Givens 
rotation  device,  a  herringbone  multiplier,  and  a  hybrid  mode  device.  Comparison  with 
experimental  data  has  also  been  included.  The  agreement  with  the  experiment  simulta¬ 
neously  verifies  both  the  electrostatic  analysis  and  the  diffraction  analysis.  Finally,  an 
efficient  configuration  for  crystals  like  lithium  niobate  has  been  identified  along  with  the 
propagating  mode  direction,  for  electrooptic  applications  involving  grating  diffraction  and 
propagating  guided-modes  in  uniaxial  anisotropic  waveguides. 
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APPENDIX  A: 


BOUNDARY  CONDITIONS  AND  FIELD  EXPRESSIONS 


The  matrix  A  of  £q.  (1)  is 

/  “ 


A  = 


«■■» 


^rntst 


e.i 

^tMwi 


—P' 


0 

0 

0 


^rnywi^aBmi 


^BBwi 

mi 

C 

^sswt 


C* 

BMW% 

^mmmi 


^ttimi 


^ammi 

13 

-  1 


0 


^aami 

^aaw  I 


(Al) 


where  »  =  e,  /,  «. 


The  boundary  conditions  of  the  anisotropic  slab  waveguide  of  Fig.  2  are: 

(a)  Film-substrate  interface  (x,  =  0) 

WfC.^wA,  (A2) 


where  the  unknown  constant  matrices  are  of  the  form  C/  =  \Fo+ ,  Fa+  *  Fo- 1  F*-]^ 
C,  =  (0,  0,  5o  _ ,  . 

(b)  Film-cover  interface  (x,  =  d) 


Wf\exp{-jk,ljd)]Cj  =WA> 


(A3) 


where  C,  =  [<70+ ,  Cs+ «  0,  0]^  and  the  matrix  C,  was  shifted  by  exp(+j^  Aed)  in  Eq.  (2) 
for  simplicity  in  the  resulting  boundary  conditions.  Combining  Eqs.  (A2)  and  (A3)  the 
following  condition  must  be  satisSed 


(A4) 


where  for  simplification  in  the  notation  [exp(-ifcoA/d)]^7**^«  =  F(j9'). 
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Using  Eq.  (2)  and  the  other  Maxwell’s  equations  it  can  be  shown  that  the  field  ex¬ 


pressions  in  the  film  region  for  a  TE  mode  are  , 

—  (■^04  ®0  4  (*■  )  "h (*•  )]  ®Xp(— *•  )»  (•^5) 

—  ^••/o['fo4  ®04  “  •^O- ®0-  )t 

Vo^tmJ  ~  ^1-^4  ®04  (*•  )  "h  ^O- ®0  -  (^»  )]  J^oP 

=  Et^f  =  fioHg^f  =  0,  (A8) 

where  ao4  (x.)  =  exp(-yA:ofc,,/oX.).  wid  Oo-(*.)  =  exp(-f-yfco*«./o*.)-  Similarly  to 
the  TE  mode  case  the  fields  in  the  film  region  for  a  TM  mode  are 

B' 

E,,j  -  {/’e  +  a£+(i.)-HF£-0£_(x,)]exp(-jfco^2.),  (A9) 

Xc««/ 

E,^f  =  (Fe  +  0£+(x,)-F,_fls-(*.)]exp{-yA;o^'x,),  (AlO) 

^oj 

rioHy,s  =\Fs^ae^{x,)-^-Fe.as.{x^)]exp{-jk;>B'z^h  (AH) 

E^^j  —  —  ^oBgrgf  —  0,  (A12) 

where  a£+  (x.)  =  exp(-j/co*;../f:X,),  and  C£_  (x.)  =  exp(+jAo  ^BW  f  E  )  . 


APPENDIX  B: 

PARAMETERS  OF  THE  GUIDANCE  CONDITION  WHEN  OPTIC  AXIS 
LIES  IN  THE  ELECTRODE  PLANE 


(Bl) 

(B2) 

(B3) 
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**  n  L3  1.4  » 


Ao£  = 


_likLfo-hj){kl,,^-b,) 


it*  -  it* 

'‘••/O  *«»/£ 


•nd  ®;p  ~  (l  P  “  (l  0^  l^tMmp^^MMmpt  ^p  ~~  ^p*mpt  Mid  p  —  Cff,S 

(corresponds  to  cover,  film  and  substrate  regions  respectively).  The  normalized  wavevector 
components  ^a««£i  k^mfot  ^d  it,./^  can  be  found  using  the  wavevector 

surface  equations.  The  ordinary  components  can  be  found  from  the  equation 


pO 


+  p  =  e,f,a. 


Because  of  inequalities  (4)  and  (5)  the  components  in  the  cover  and  the  substrate  are  pure 
imaginary  and  their  sign  should  be  chosen  in  order  to  have  exponentially  decaying  fields 
in  these  regions.  The  extraordinary  components  that  appear  in  the  film  and  the  substrate 
regions  are  given  by  the  extraordinary  wavevector  surface  equation 

+^'^)  +  ("£p P=/.  «.  (B7) 

where  e,.  is  the  direction  cosine  of  the  optic  axis  along  the  waveguide  a,  coordinate 
axis.  The  substrate  component,  because  of  inequality  (4),  b  pure  imaginary  and  special 
attention  should  be  taken  in  the  choice  of  its  sign  in  order  to  have  exponentially  decaying 
fields  in  the  substrate. 

The  parameters  where  p,  q  =  E  Bit 

.  _  .  - 1  /  ~  ~  “I  #pg\ 

1  Ao -»,))(.,/„)  *..,o  /  '  ' 

V  ■  A,  +  lAc,/(lt2...  -  ».))(»./.,)  k..,,  I 

where  q  —  O,  £,  and  the  rest  of  the  parameters  have  been  previously  defined.  The  param¬ 
eters  Go ,  Ge ,  Ho ,  Hs  ^  ffo ,  9e  *  ho ,  hg ,  ao ,  and  Og  are,  G,  =  (l/2)ll  +  *2,,,/|fc.«e|’l^'’, 

p,  =  tan-‘(*.././lfc.-.l).  B,  =  +a*nS./|k...|*l>/»/2|A:2.„  -6,1,  h,  = 
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tan“*(a/n2,^/c;ifc., and  a,  =  where  q  —  E  and  d  is  the  film 

thickness. 


The  expressions  for  the  fields  in  the  film  region  are 


^-/=  E  ^-^[F,^fl,.,(x.)-F,.a,.(z.)]exp(-y*o/?'*.),  (BlO) 

E. .,=  T[F,^a,^(x.)  +  F,.a,.{x.)\ap{-ik^ff'z.),  (Bll) 

f»o.a 

F. ,/ =  E  +  ^«- «.-(*- )]<«p(-yko;9'i.),  (B12) 

,.0.£ 

rioff.m/  E  ®*+  exp[-jkofi'z, ),  (Bl3) 


t«0.£ 


VoH^wj  =  E  )  -  -^-  o«-  )]  exp(-jfeo)g'g» ).  (B14) 

VoH.m/ -  E  (®15) 


1*0, s 


The  expressions  for  the  Bot  Be,  and  that  appear  in  the  expressions  for  Fb  + ,  Fo - , 
Ff  4 ,  and  F^  _  in  the  case  of  the  optic  axis  lying  in  the  plane  of  the  electrodes  are 


Bo  ~  ((Foo  cos  d>oo  —  Foe  Q cos  4*ob 

(Foo  sin^oo  "■  FoeQsin^oa)*)^^** 

Be  =  ((Fjb  o  cos  o  ~  Eg  eQ  cos  £  )*  + 

(Fbo  *in^£o  *“  F££Q8m^££)’]*^’, 
.  _  -1  f  Fqo  sin^oo  ~  FqeQsir^oe  \ 

^  '  Foo  C08 ^oo  —  Fo£  Qcos  ^o*  ^ 

«  -I  f  F£o  8in^£o  —  F££^sin^££  \ 

*  t  Eg  o  cos  ^£  o  *“  Eg  £  Q  cos  ^g  g  J 


(B16) 


(B17) 

(B18) 

(B19) 


?£.  Fqo  8in(ao  -  hp  +  ^qq)  Hg Ego  Bin(a£  —  hg  4-  ^£o) 
BqEoe  8in(oo  —  fcp  +  4os)  +  BgEgg  6in(a£  —  hg  ^be) 


(B20) 
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APPENDIX  C: 


POYNTING  VECTOR  COMPONENTS 


1.  Optic  Axis  Perpendicular  to  the  Electrode  Surface 

The  i  component  of  the  Poynting  vector  for  this  optic  axis  orientation  is 

=  ~ — +  cos(2fcofc,,/o®»i  +  2^/,o)]l/iO+  |*P*l^o- 
Vo 

for  the  T£  modes,  and 

S,wif  —  ~ — ^(1  +  cos(2fcoA:,«/£X,i  +  2^/,£)]|/i£+ I’, 

Vo^MMmf 


for  the  TM  modes  [where  rjo  = 


2.  Optic  Axis  in  the  Plane  of  the  Electrode  Surface 


For  this  case  the  component  of  the  Poynting  vector  is 

S.../  =  f  {BS  (1  +  +  Bi  (1  +  -^)+ 

c* 

B%{l - cos(2fco^,«/o**i  +2do)+ 

Bl{l  -  cos(2fcofc,.,*x,i  +2^^)+ 

^•twf 

4Ro^£[cos(*o<:.«/o*-<  +  «?o)cos(*o*../s*.i  +  «>x)+ 
CoCk  Bin(fcofc../o*«<  +^o)«»n(A:ofc.,/*x,i  +  tf*)]} 
p*i/io+ n^o-  I’l 


where  C,  =  Kw ”  */)»  9  = 
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FIGURE  CAPTIONS 


1.  (a)  Three-dimensional  geometry  of  the  diffraction  of  guided-waves  in  uniaxial  an¬ 
isotropic  slab  waveguides  by  interdigitated-electrode  induced  phase  gratings,  (b)  A 
top-view  of  the  same  geometry  that  shows  the  various  diffraction  directions  of  the  or¬ 
dinary  and  extraordinary  components  of  an  incident  mode,  for  each  diffracted  order. 
In  this  figure  only  the  0-  and  +  1-orders  are  shown. 

2.  The  decomposition  of  a  guided  mode  into  four  plane  wave  components.  0-1-  and 
O—  are  the  two  ordinary  plane  wave  components  while  E-\-  and  E—  are  the  two 
extraordinary  plane  wave  components. 

3.  Decomposition  of  the  guided  mode  diffraction  problem  into  four  grating  diffraction 
subproblems. 

4.  Planar  slices  through  wavevector  surfaces  for  the  case  of  the  optic  axis  being  per¬ 

pendicular  to  the  electrode  surface,  (a)  zz^  plane  intersection  at  y.  =  0,  (b)  xy 
plane  intersection  at  z  =  0,  (c)  xy  plane  intersection  at  x,  =  »  correspond¬ 

ing  to  the  ordinsiry  waves  0-1-  and  O—  (TE  modes),  and  (d)  xy  plane  intersection 
at  X,  =  ±kg^jB,  corresponding  to  the  two  extraordinary  waves  E+  and  E—  (TM 
modes). 

5.  Top-view  (xy  plane)  of  the  diffraction  geometry  in  the  case  of  the  optic  axis  perpendic¬ 
ular  to  the  electrode  surface.  The  0-order  and  the  t-order  forwzLrd-diffracted  directions 
are  shown,  (a)  TE  incident  mode  (O-f,  O—  waves),  and  (b)  TM  incident  mode  {E+, 
E-  waves).  The  directions  shown  are  valid  for  a  negative  birefringent  material.  In 
the  case  of  positive  birefringent  material  the  relative  ordering  of  the  diffraction  angles 
is  reversed. 

6.  PlanaLT  slices  through  wavevector  surfaces  for  the  case  of  the  optic  axis  being  in  the 
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electrode  plane,  (a)  zz„  plane  intersection  at  y.  =  0,  (b)  zy  plane  intersection  at 
z  =  0,  (c)  zy  plane  intersection  at  z„  =  corresponding  to  the  two  ordinary 

plane  wave  components,  and  (d)  zy  plane  intersection  at  z«,  =  ,  corresponding 

to  the  two  extraordinary  plane  wave  components. 

7.  Top-view  (zy  plane)  of  the  diffracted  geometry  in  the  case  of  the  optic  axis  lying  in  the 
electrode  plane.  Only  the  0-order  and  the  i-order  forward-diffracted  waves  are  shown 
for  (a)  general  incidence  and  optic  axis  orientation,  and  (b)  for  Bragg  incidence  and 
optic  axis  orientation  restricted  along  the  z  or  the  y  axis. 

8.  zy  plane  intersections  of  the  wavevector  surfaces  at  z„  =  that  correspond 

to  the  ordinary  waves,  and  at  z„  =  that  corresponds  to  the  extraordinary 

waves.  The  Bragg  condition  is  also  satisfied,  (a)  Optic  axis  along  the  z  axis,  and  (b) 
optic  axis  along  the  y  axis. 

9.  The  general  relation  of  the  complex  coefficients  fio+ ,  /,o-  *  /<s+ »  and  /<£_  ,  in  the 
complex  plane.  Coefficients  ratios  with  respect  to  /,©+  are  shown. 

10.  zy  plane  intersections  of  the  wavevector  surfaces  for  the  case  of  the  optic  axis  being 
perpendicular  to  the  electrode  surface,  and  for  a  slanted  electrode  configuration,  (a) 
A  TE  incident  mode  (ordinary  waves),  and  (b)  a  TM  incident  mode  (extraordinary 
waves). 

11.  zy  plane  intersections  of  the  wavevector  surfaces  at  z,  =  ±fc,,/o  (ordinary  waves) 
and  z,  =  ±/k,,/E  (extraordinary  waves)  for  Bragg  incidence  and  (a)  optic  axis  ori¬ 
entation  parallel  to  the  grating  vector,  and  (b)  optic  axis  orientation  perpendicular  to 
the  grating  vector. 

12.  Top-view  (zy  plane)  geometry  of  a  Givens  rotation  device.’’ 

13.  The  diffracted  electric  fields  of  the  Givens  rotation  device  as  a  function  of  the  applied 


2 


electrode  voltage.  The  dotted  lines  correspond  to  rigorously  calculated  results  while 
the  continuous  lines  correspond  to  a  best  cosine  or  sine  fit.  (a)  0-order  forward- 
diffracted  fields,  and  (b)  -I- 1-order  (Bragg)  forward-diffracted  fields. 

14.  (a)  The  top-view  of  the  herringbone  multiplier,  (b)  The  wavevector  diagram  (xy  plane 
intersection)  that  corresponds  to  the  herringbone  multiplier.  The  incident  wavevector, 
ki ,  the  Bragg  diffracted  wavevectors,  and  S2 ,  and  the  grating  vectors,  and  K2 
are  shown. 

15.  The  efficiency  of  the  —2-order  forward-diffracted  mode,  the  predicted  efficiency  from 
the  approximate  formula  DE^r  =  K,ViV2,  and  the  percentage  error  [DE^r  -  DE)  x 
100/Df?,  as  functions  of  the  applied  voltages  (V  =  Vi  =  Vj). 

16.  Diffraction  efficiencies  of  the  0+,  O— ,  and  E—  forward-diffracted  waves  as  a 

function  of  the  applied  voltage  for  (a)  the  0-order,  and  (b)  the  Bragg-order. 

17.  Diffracted  phases  of  the  0+,  O— ,  JS-f ,  and  E—  forward-diffracted  waves  as  a  function 
of  the  applied  voltage  for  (a)  the  O-order,  and  (b)  the  Bragg-order. 

18.  Diffraction  efficiencies  of  the  O-order  forward-diffracted  mode  (dashed  line),  and  of  the 
Bragg-order  diffracted  mode  (continuous  line)  as  a  function  of  the  applied  voltage.  In 
these  plots,  the  computed  solutions  are  given  by  the  solid  dots.  The  continuous  and 
the  dashed  lines,  calculated  by  cubic  spline  interpolation,  have  been  provided  for  ease 
of  interpretation. 

19.  The  theoretical  (continuous  line)  and  the  experimental  (circles)  mode  diffraction  effi¬ 
ciencies  of  the  interdigitated-electrode  correlator*'  as  a  function  of  the  applied  voltage. 
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